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' Abstract 



We present a framework for the analytic calculations of the hierarchical wave 



o 

CJ ■ functions and the composite fermion wave functions in the fractional quantum 



■ Hall effect on the sphere by using projective coordinates. Then we calculate 

>< 

■ the overlaps between these two wave functions at various fillings and small 
numbers of electrons. We find that the overlaps are all most equal to one. 
This gives a further evidence that two theories of the fractional quantum Hall 
effect, the hierarchical theory and the composite fermion theory, are physically 
equivalent. 

PACS: 73.40.Hm, 73.20.Dx,03.65.-w,03.80.+r 



1 



Typeset using REVTgX 



I. INTRODUCTION 



The fractional quantum Hall effect (FQHE) at the Landau-level (LL) filling fraction 
V = 1/m with m an old integer is very well described by Laughlin's theory.0l The Laughlin 
wave function is a very good approximation of the exact ground state of the quantum Hall 
effect (QHE) at v = 1/m. However for the FQHE at u ^ l/m, there exist two well-known 
theories (notice that we will only consider the case that the electron spins are polarized in 
this paper). One is the hierarchical theory. The states at z/ 7^ l/m are formed due to the 
condensation of the anyonic quasiparticles of Laughlin states.iHlSl The trial wave functions 
constructed from this theory are called as the hierarchical wave functions. Another theory 
is based on the composite fermion (CP) approach proposed by Jain,0 where the FQHE is 
due to the integer QHE of the composite fermions (CFs) (electrons bounded with an even 
magnetic flux quanta). The trial wave functions constructed from the CF theory are called 
as the CF wave functions (or Jain's wave functions) The overlaps of the exact states with 
the hierarchical wave functions and the CF wave functions are both excellent. It has also 
been shown that two theories predict the same topological excitations at the same u .Billi 
The two theories must be physically equivalent if they both describe correctly the physics 
of the FQHE. Thus it would be very interesting to study the difference and equivalence of 
the two theories. 

In this paper, we present a framework for the analytic calculations of the two wave 
functions on the sphere by using projective (or stereographic) coordinates on the sphere. 
There are several advantages of using spherical geometry. As it is a compact surface, there 
will be no edge state to be worried if we are only interested in the bulk state. Also the system 
has rotational invariance symmetries. On the torus, though the system has translational 
invariance and no boundaries, the hierarchical wave functions are very difficult to calculate 
and quite complicated due to its nontrivial topology,^! and we do not even know how to 
construct the CF wave functions with the correct center coordinate degeneracy on a torus. 

Because the states considered in the FQHE are restricted to the lowest Landau level 
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(LLL), the wave functions are only dependent on holomorphic coordinates (polynomials of 
the holomorphic coordinates) on the sphere. Therefore it is possible to use only holomorphic 
coordinates to do all calculations. To compare the two types of hierarchical wave functions 
is the same to compare the two polynomials of holomorphic coordinates on the sphere. We 
note that our ultimate goal is to expand those wave functions in polynomials and calculate 
the overlaps of two wave functions or physical quantities (for example, the density- density 
correlations) at an arbitrary number of electrons by the method (Jack polynomials method) 
used in studying the Calogero model.lii We do not know how to do it at the moment, and 
further progresses on it will enhance our understandings of the theories of the FQHE. 

We organize the paper as follows; first we review the Landau level problem on the sphere. 
A self-contained derivation of eigenstates of an electron on a sphere with a monopole field 
by using a simple geometric argument and projective coordinates is given in the appendix. 
Then we show how to classify the many-body eigenstates of the angular momentum in the 
LLL. We then construct the wave functions based the theory of the hierarchical states and 
the theory based on the CF picture. The wave functions constructed in this paper are easy 
to handle in the practical calculation. Finally we calculate the overlaps of the hierarchical 
wave functions and the CF wave functions at various fillings v and some small numbers of 
electrons. 



II. QUANTUM MECHANICS ON THE SPHERE 

The electrons are constrained to move on the surface of a sphere of radius R having a 
magnetic monopole in its center. The total magnetic flux AnR^B must be an integer multi- 
ple = 25 of the magnetic flux quantum 0o = ^hirc/e according to the Dirac quantization 
condition. Therefore, the sphere radius R is equal to S^^'^Iq, where Iq = ij^Y^"^ is the mag- 
netic length. The eigenstates of an electron are given by monopole spherical harmonics.iJii 
First, we briefly review the old method to derive the wave functions of the Landau levels 
(LLs), then rederive them by using Algebraic Geometry. 
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For simplicity, we take units h and c equal to one in the following formulas. The Hamil- 
tonian of a single electron of mass rrie is given by if = 2^(-P* + eA)^. However, since the 
electron is confined on the spherical surface, one shows, 

where A = r x (P + eA), Uc is the cyclotron frequency, P = — ^V , V x A = Bfl, and 
n = r/R. 

The components of A obey the commutation relations [Aj, Aj] = ieijk{Ak — SQk)- The 
angular momentum operators L = A + Sfl, and their commutation relations are [Li, Lj] = 
ieijkLk- Since A is normal to the surface, we have f2 ■ A = A • = 0, and L ■ f2 = f2 ■ L = 5*. 
Using these equations, the relation |Ap = |Lp — can be obtained. Thus the eigenvalues 
of |Ap can be deduced from the usual angular momentum algebra |Ap = |Lp — S*^ = 
L{L + 1) — S"^, L = S + n, n = 0, 1,2..., and the eigenstates of the Hamiltonian are the 
eigenstates of |Lp and L3, and they are given by monopole spherical harmonics. We choose 
a gauge field A = —-^^^-^^'^1 of which the singularity lies on the north pole (we choose 
a different gauge from the one used in Ref. |^). The wave functions at the LLL are given by 

^S+m^S-m^ (2) 

where m = —S, —S + 1, - ■ ■ , S, and 

u = cos(^^)e^^, V = sin(J^). (3) 

All wave functions of the LLs can be derived by this way,lii and we will not repeat this 
derivation here. In the following, all eigenstates will be obtained by using projective 
coordinates,lli0 and the method developed in Ref. [18 . 



The projective coordinates are given by 2; = 2_Rcot le*"^ and its complex conjugate z. We 
will take R = 1/2 for simplicity. The measure on the sphere is / jff^^- The Hamiltonian 
of Eq. (|l|) in projective coordinates is now written by the following formula,0 

H = —(1 + zz)\P, + eA^m + eA,), (4) 
rrip 



where 



OZ oz 2 1 + zz 



(5) 



and (f) is the flux (in the unit of the fundamental flux ^o) out of the surface. Note that the 
Hamiltonian given by Eq. (H) (we call this Hamiltonian as H' in the appendix) is different 
from the one given by Eq. (|l]) by a constant. 

The ground states can be determined from the solutions of the equation {Pz + eAz)'ip = 0, 
and they are (unnormalized) 



where / = 0, ■ ■ 
the appendix) 



(6) 



1 + zz)'2 

At any Landau levels, the eigenstates (unnormalized) are given by (from 



^n,i = [d, + (f + ^)djng][dz + (f + 2)djng] 



B 



(0) 



x[dz + {- + n - l)djng]^':> 



(7) 



where 



1 



[i + zzy 



B = 0/2, 



(0) 



1, z, . . . , z , . . . , z 



(l>+2n 



(8a) 
(8b) 
(8c) 
(8d) 



Under any finite rotations, z coordinate is transformed as z' = The rotation matrix 



R 



yC d j 



is generated by the rotations along the three Cartesian axes. 



1 

V2 



V2 



' (l + cosa)i/2 z(l -cosa)i/2^ 
^ i(l - cosa)^/^ (1 + cosa)^/^^ 

' (l+COS/?)l/2 (l-COS/?)V2^ 
^-(l-COS/5)l/2 (l + C0S/?)l/2^ 



(9a) 



(9b) 



Rz — 



exp(i7/2) 

exp{—ij/2) 



\ 



J 



(9c) 



The rotational invariance of Hamiltonian is shown by the identity: 

OH{z')0-^ = H{z), 

where 

= (^_)l. 
The wave function is transformed under rotations as 



(10) 



(11) 



' cz + b 

We hst some useful relations when we do a finite rotation on a many-body wave function 



(12) 



d(^Zi, Zj 



"I ~j 



VI + ZiZi^l + ZjZj' 



Zi — Zj, and 1 + ZiZi are transformed under the finite rotation as 



Zi Zn 



{czi + d) {czj + d) ' 

1 + ZiZi 



[cZi + d) {cz + d) 
Finally, the angular momentum operators for N electrons are 

N 
i=l 

1 ^ d 

N 

Jy — Jy{^) 



(13) 

(14a) 
(14b) 
(14c) 



(15a) 



i=l 



N 



2E[(i + -n|: + (i + ^D^ + t(-.-^^. 



dz 2 



(15b) 



TV 



i=l 
N Q 



(15c) 
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III. PROJECTIONS AND ANGULAR MOMENTUMS IN THE LLL 



The FQH state is restricted to the LLL. In this section, we will discuss briefly how to 
project states to the LLL on the sphere (see Ref. |T^, and for the detailed discussions in the 
case of a plane or a disk, see Ref. ^), and how to find the eigenstates of angular momentums 
when the particles are restricted to the LLL. Note that the construction of the CF wave 
functions involves the higher LLs, we need to project the wave functions to the LLL (see 
Sec. 0). 

The normalized states with flux in the LLL are 

and / = 0, 1, 2, ■ ■ ■ , 0. The projection operator to the LLL is P = X)/ 1^ >< Ai ^^'^ can be 
written also in the following form, 

Pi,{z, -z)= f ^^^Giz, w)i^{u, CO), (17a) 

n( \ 0+1 {\^zwf 

""^ = 2vr (l + z-zy/^il + wwr/^ - ^'^^^ 

For the many-body wave functions, P (or G) is equal to Dili Pi (or Hili Gi) where Pi is the 

projection operator of the i-th particle and is the number of particles. 

If the state is not in the LLL, the anti-holomorphic coordinate z will appear. Typically, 

it appears as 

^ = (TT^Iw^' 

and Pip is equal to 

(0 + l)!(/ + O!(0 + j -/-2)! z' 



(19) 



/!(0-/)!(0 + j + l)! (l + zz)^/^' 
On the sphere, if the interactions between electrons are rotationally invariant, the eigenstates 
of the many-body Hamiltonian should be also the eigenstates of rotational operators 
and Jz- The FQH ground states on the sphere are rotationally invariant and are non- 
degenerated. In order to find the ground states, we can thus use the rotational invariant 



states to diagonalize the Hamiltonian. As the number of all possible rotational invariant 
states is much less than the number of all possible states, it is thus much easier to find the 
ground states by using the rotational invariant states to diagonalize the Hamiltonian than 
by using all possible states. It could be also interesting to find the eigenstates of 7^ 
(which are not rotationally invariant). The excited states in the FQH are not rotationally 
invariant. For Fermi-liquid-like systems in a half-filled Landau level, one can have ground 
states which are not rotationally invariant .0 

Now we are going to find the many-body wave functions on the LLL which are the 
eigenstates of and Jz- In the LLL, the many-body wave functions \I' have the form 

N ^ 

^ = n T^Pi^i^ ^2, ■ ■ ■ , ^7v), (20) 

i=l (1 + ZiZi)2 

where 2:2, ■ ■ ■ , -Zat) is an anti-symmetric holomorphic function. When Jj^ = + 

iJy, J_ = Jx — iJy, Jz act on we have 

TV ^ N g 

i=l (1 + ZiZi) 2 i=i oz^ 

i=l (1 + ZiZi) 2 i=i OZi 

^1 ^ d N6 

■'•* = n7--^i(i:^.^)-^if- (2i<=) 

i=l (1 + ZiZi) 2 i=i OZi / 

Thus the projected J operators are: 

N p) 

N fj 

J'-, = T.-^'g^ + <P^^^ (22b) 
i=l ^ 

^: = (E^.S-)-^. (220) 

1=1 * 

where they act only on F. The angular momentum eigenstates of the many-body wave 
functions restricted to the LLL can be obtained by solving 

J'_F{-J) = 0, (23a) 
J'zF{-J) = -JF(-J), (23b) 
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where F{—J) is the lowest weight eigenstate with weight — J. Other states can be obtained 
by applying repeatedly on F{—J)'. Eq. (|23| ) leads to 

J:^^-^H-J) = - JM-J), (24a) 

j:^Fi-J)=0. (24b) 

i=l * 

The first equation in Eq. (^) means that F{—J) is a homogeneous polynomial with degree 
^ — J. As F{—J) is an anti-symmetric function of holomorphic coordinates, it can be 
factorized as F(— J) = Y[ii:jizi — Zj)F'{—J). One can check that 

/ F(-J) = n(-2i - z,)J'_F'{~J), (25a) 
J:F(-J) = ^ , ^ F(-J) + n(^. - ^.)^:i^'(-^). (25b) 

Thus F'(—J) is a symmetric function with degree L = ^ — J — Mi^nl)^ the power of 
every coordinate in F'{—J) shall be less or equal than 0' where </)' = (/) — {N — 1). By using 
Eq. ( p4|) and Eq. (^), one finds that F\—J) satisfies the conditions: 

J'_F\-J) = 0, (26a) 
J'.F'{-J) = (^ - J)F'{-J). (26b) 

Define symmetric polynomials af. 

N N 

p{z.) = m^-^^)=j:i-^y^^^''-\ (27) 

i=l 1=0 

where 

AT N 

(To = 1, cTi = ^ Zi, ■ ■ ■ , (Tat = Jl Zi. (28) 

1=1 i=l 

F' can be expanded as 

N 

T.C{s^)U^r, (29) 

Si i=l 

where Sj is a non- negative integer. By using Eq. (^), we get equations which C{si) and Sj 
must obey. One of them is 



The condition 

Es.<0' (31) 
1=1 

must be satisfied in order that the wave function is normalizable. C(sj) shall also satisfy 
the equation 

J-F'i~J) = T.C'isdU^:' = ^^ (32) 

s' i=l 

where C'{s[) is a linear combinations of C(sj), and it shall be equal to 0. Because [J^ , J_] = 
— J_, s[ shall satisfy 

TV 

Ez.: = L-l. (33) 

i=l 

Apparently Y^^=i s'i shall be also less or equal than cp' . The number of linear independent 
solutions for C(sj) is equal to the number of solutions of Eq. (^) minus the number of 
solutions of Eq. (^), and it is also equal to M( J, N, 0'), which is the number of states with 
spin J. 

The generation function for the number of solutions of Eq. ( pOj ) or Eq. ( P5D is 

G(*..) = nT^. (34) 



1=1 



The number of solutions of Eq. ( pO] ) is equal to the sum of the coefficient of term Pq^ with 
< J < 0' in G{t, q). Thus M(J, A^, 0') is equal to 

^^G(t,g)(l-^)f:i. (35) 

We can also use a generation function of one variable,iil 

Git) = ^ ^^'=' . (36) 

The number of solutions of Eq. ( POD is then given by the coefficient of of function G{t). 
Thus 
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The asymptotic behavior of M( J, A^, 0') can be obtained by using the steepest descent 
method.ii When L and 0' are both large, M{J, N,(f)') is equal to exp{S{J, N,(f)')) approxi- 
mately, and exp{S{J, N, 0')) is determined by the following equations, 

^ exp{-p) ^ 1 ^ r{N+4,')p ^ rNp 

1 — exp{—p) Jo Jo 

+ f%u ^'^P(-P^) (38a) 
JO 1 — exp[—pu) 

SU N, 0') = Lp + ln(l - exp{-p)) + -( / - / 

p Jo Jo 

+ r ^)du ln[l - exp{~u)] . (38b) 

JO 

We list the number of rotational various fillings in Table |. 

IV. HIERARCHICAL WAVE FUNCTIONS 

In this section, we will discuss the construction of the hierarchical wave functions. 
The quasiparticles satisfy fractional statistics, and the condensation of quasiparticles gives 
rise to the FQH state with z/ 7^ 1/m. 

Define 

N 

= Ui^iVj - UjVir, (39) 
i<j 

where m is a positive integer. For = ^ with m being an old positive integer, the FQH wave 
function will be \E'm (the Laughlin wave function). The flux is equal to 0^ = m(N — l).i 
Or in the projective coordinates, it is 

N 

^m = Y[dizi^^jr- (40) 

i<j 

The Laughlin wave function with the presence of quasiparticle excitations is given by acting 
the quasiparticle excitation operators on the original Laughlin wave function. The quasi- 
particle excitation operator is given by 
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N 

A\a^l3) = Y\_iP'^i ~ Q^"^*)) (quasihole) (41a) 

i=l 

^ - d _ d 

Aia, [3) = W{[3— a——), (quasielectron) (41b) 

f^i oui dvi 

where a = cos |e*'^, and P = sin | are the quasiparticle coordinates. In the projective 
coordinates, the operators of the quasihole excitation and the quasielectron excitation are 
given in the following form, 

N 

A^icu, io)'^m{zi) = n dizi, Uj)<iJm{Zi), (42a) 

1 1 



N 



X 



Yliil + ZiUj)d,^ - (l)mio]F^, (42b) 

i=l 

where a; ,ci; is the projective coordinates of the quasiparticle, and Fm{zi) = Y[ii:j{zi — zj)"^. 
The flux in the presence of a quasielectron (quasihole) is 0^ — 1 (0^ + 1)- 

The slightly entangled appearance of A^u) hides, indeed, a form which is analogous to 
A^{uj). To unveil the similarities, one can show that 

TV 

Pi^^-l,Zi)Y[d{Zi,Cu)^m (43) 
i=l 

gives the wave function of the Laughlin state in the presence of a quasihole as that in Eq. 
(^). P{(f),Zi) (here = (pm — 1) projects the wave function to the LLL with flux with 
respect to coordinates Zi. Thus the construction of the hierarchical wave functions due to 
the condensation of quasielectrons will naturally involve higher Landau levels as in the case 
of the CF wave functions (see the next section). 

Instead using A{uj , u)"^ m{zi) , we can also create a quasielectron excitation using 
^rn-2AD^E where D = [^i(^i)]^. AD is here equal to P(02 - l)Uf=i d{zi,u)D. We 
call (\E'i)™~^y4D as a wave function by the hard core construction. 

In the case of many quasiparticle excitations, the operators of excitations are 

= Il^^i^J^^j)^ (44a) 



"9 
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= n^(^j'^j-)- (44b) 



When An^ acts on one can show 



xA'F^izi), (45a) 



where 



= n ^v) 



N 
i=l 

N 
i=l 

-{(j)m - Nq + 2)uNq-l] 



N 

X + Ziu;j)d^^ - {(f)m -j + 

i=l 



N 

X ]J[(1 + ZjLJi)d,^ - 0m^i]. (46a) 

i=l 

One shall be careful about the ordering of A'{j) in Eq. (^6|). A'(j) in Eq. (^Bf ) is defined as 

N 

^'(j) = n[(l + -^^^jO^.. - (0™ - J + (47) 

j=l 

As in the case of a single quasihole excitation, the wave function in the presence of Nq 
quasiholes can be also written as 

N Nq 

P{(p)l[I[d{z^,COa)'fm, (48) 

i=l a=l 

where = $m — Ng. 

To construct the hierarchical wave functions, we shall normalize the Laughlin wave func- 
tions in the presence of quasiparticles. One can show that \E'c,q = [^i(^^Q)]^^"^^]v,^m5 ^e,q = 
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[^i{uJa)Y^"^AN^'^m, or for the hard core constructed wave function, i{uJa)Y^"^'^ m-2AD , 

are normahzed.ii The hierarchical states are obtained if the quasiparticles are also condensed 

to the Laughlin states. The wave function of quasiholes is = ['^i{iUaY'^~^"^ , the wave func- 

_ j_ 

tion of quasielectrons is \E'g = [\E'i(ti;a)]^^ and p is a positive even integer. Quasiparticles 
satisfy fractional statistics,! and the wave functions here are in singular gauge which shows 
fractional statistics explicitly. The hierarchical wave function for the electrons is then given 
by the following formula, 

\ (1 + UJatOaf 



or explicitly, 



N Ng 



Mm,p) = / n ,,7 n n ^(^^,^«) 

a=l{^ + i=la=l 

x\^,M\^/^^p,{u^) (50) 

is the hierarchical wave function due to the condensation of quasiholes, and the hierarchical 
wave functions due to the condensation of quasielectrons by the non-hard core construction 
and the hard core construction are given in the following formulas. 



N Nq 

^ n n d{zi,Ua)'^p^{uJa), (51a) 



i=l a=l 



^'e,hard(m, -p) = ^™_2(^i)^(02 - Ng, Zi) J Y[ 



a=l (1 + ^a^a)^ 

N Nq 

x^2(^*)n n d{z^,tOe,)'^p,M. (51b) 

i=l a=l 

We also require the wave functions above being rotationally invariant. This requirement 
leads to 

m{N - 1) + ^2Ng = 0, (52a) 
P2{Ng-l) = N. (52b) 
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^2 = ±1 in the case of the condensation of quasiholes and quasielectrons respectively. And 
the Landau level filling fraction v is equal to 

(53) 

For m = 1 and ^2 = 1, the filling z/ = — ^- is equal to the filling of the charge conjugate 

m+ — 

P2 

state, 1 — Actually, the wave function \l/e(m,p) is also the charge conjugate of the 

Laughlin wave function at filling u = and this shows that the construction of the wave 
function is consistent with physical picture. When m 7^ 1, we notice that, in the formula 
for \E'e(m,p), we can not do the integration exactly due to the term |\I^i(ci;q) p/™. We can 
approximate the trial wave function \I^e(^5P) by omitting |\E'i(u;Q,)p/'", and it becomes 

J [i- + UJaUJa) 

N Nq 

X J] Jl (54) 

i=l a=l 

The wave function written in Eq. ( [5^ ) is still rotationally invariant and we are able 
to integrate it. When m = 1, the formula for \I'e(m,p) in Eq. (^) is integrable. When 
m = 1, we find that the overlap between the wave functions given in Eq. ( pD] ) and Eq. 
is excellent for a small number of electrons. In Ref. ^ it was also found that the overlapping 
of the wave functions given by Eq. (|54D with the exact ground state of the FQH is all most 
equal to one for a small number of electrons. We note that the wave functions calculated in 
Sec. are based on the formulas written in Eq. (0) and Eq. (|5^). 

In the formula for \l/e(^, —p) or \l/e,hard('^; ~p), we note that one can do the integration 
first, and then the projection, or vice versa. In Sec. 0, the overlap between these two wave 
functions will be calculated for a small number of electrons and it is found that the overlap 
is all most equal to one. 

We will call the above hierarchical states as the 2nd-\evel hierarchical states, and the 
Laughlin states as the Ist-level hierarchical states. The higher-level hierarchical states 
can be built in the similar way.BiSi We denote the k'th-\eve\ hierarchical states by 
{Pi,^2P2,(,3P3, ■ ■ ■ ,(,kPk), where pi is an old positive integer, pi,i ^ 1 are even positive 
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integers, and = ± indicate the quasihole condensation and quasielectron condensation 
from parent states. 

For the higher-level hierarchical wave functions involving the condensation of quasielec- 
trons, we can make a further simplification. We take {pi, —p2, 
— ps) as an example. The wave function for this state is 

a=l ' ^Q^a) i=la=l 
J j-Ji (1 + <<)2 

N2 

X UUdi^l^D'^Psi^l)^ (55) 

a=l (3=1 

where Ni is the number of electrons, N2 is the number of quasielectrons of the Laughlin 
state (pi), N3 is the number of quasielectrons of the hierarchical state (^1,-^2), 
a;^ are the coordinates of quasiparticles of the two types respectively. We can prove that 
P{Ni,u)l^) can be drooped inside the formula. Thus the wave function can be written as 



= P(0, / n n n t"i n2ri ?t^^^2 
J „=i (1 + ^M)^ (1 + ^j^W 

Ni N2 N2 N3 

xY[l[d{z,,ol)l[l[d{ulul). (56) 

i=l a=l a=l 13=1 

The wave function in Eq. (|56D is quite similar to the wave function constructed in Ref. |^, 



e 



Ni N2 1 N2 N3 -r 

><UU^r^^UU^r^y (57) 

i=l a=l "l^i' ^a) a=l (3=1 "l^a' ^(3) 

However it is difficult to handle Eq. (|57D in the practical calculation due to the singularities. 
Finally by requiring the rotational invariance of the wave function (p6D or (0), one gets 

Pl(iVi-l)-iV2 = 0, 

N,-p2iN2-l) + N; = 0, 

N2-p3iNs-l)=0, (58) 
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and Eq. (^) implies that the filhng of the FQH state is equal to 

1 



(59) 



P3 

We point out that the wave function proposed in Ref. |^ had been also constructed on the 
torus.lll It would be very interesting if we can generalize the construction of the wave function 
to the torus. 

V. COMPOSITE FERMION WAVE FUNCTIONS 

The CF theory of the FQHE has significantly advanced the understanding of the FQHE 
recently.0 The FQHE is due to the integer QHE of the CFs, where a CF is the bound 
state of an electron and an even number of vortices. We will discuss in this section how to 
calculate the CF wave functions in our framework. 

Jain proposed that all trial wave functions of the FQHE (note again in this paper the 
spin is polarized) can be obtained by using two operations, D and C, respectively composite 
fermionization and charge conjugation, on the wave functions of the integer QHE of the CFs. 
For example, the trial wave function of electrons at u = n/{2n + 1) can be written as PDxn, 
where Xn is the wave function of the CFs which fill completely the first n Landau levels with 
flux (j)* {P is the projection operator to the LLL as in the previous sections). The flux of the 
state PDxn is equal to 2(iV-l) + 0* where (p* = ^-n. We can also use ^iP(0-iV+l)^iXn 
as the trial wave function and we call this wave function as the wave function by the hard 
core construction. The charge conjugation of PDxn (or \I'iP(0 — + l)\l/ix„) is then trial 
wave function at u = 1 — tttt = tt-tt- The trial wave function at other fillings can be 

271+1 271+1 ° 

obtained by acting repeatedly D and C on PDxn {^iP{'P — N + l)\E'ix„) (each state can 
be obtained only in a unique way in this picture except the ordering of operator P) . 

Xn is given by the determinant Xn = det{ilJs,k{zi)), where s = 0, l,---,n — 1, k = 
(p* + 2s + 1, i = 1,2, ■ ■ ■ , N, N = ncj)* + n^. det{ips,k{zi)) can be simplified and it is given by 
the following formula. 
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Xn 



1 



Zl 



Zl 

ZlZi 



1 

Z2 
Zl 



N'-l 
Z2 



Zl 

Z2Z2 



Z\Zi Z2Z2 



Zl 



Zl Zl 



Z2 



Z2 Z2 



-n-\N'-\ n-lN'-l 
Zl Zl Z-i Z2 



1 

Zn 

zn 



N'-l 

Zn 



Zn 

ZnZn 



iz ^N'-l 
ZnZ^ 



Zn 



Zj^ Zn 



n-l N'-l 

Zn zpf 



n 



, (60) 

■=i {1 + ZjZj)^+--^ 

where N' — N/n — (f)* + n. We divide N electrons into n groups. The set of the origi- 
nal coordinates Zi can be mapped to Zs,k with s — 0, 1, ■ ■ ■ , n — 1, k = 1, 2, ■ ■ ■ , N'. The 
determinant is proportional to 



n-l 



Xr 



AN n [ef 



s=0 



X 



AT 

n 



1 (1 + ZiZi) 2 



where 



N' 



k=l 

N' 

— H d{Zs^ki — Zs,k2)l 

fel<fc2 



(61) 



(62a) 



(62b) 



and AN is the anti-symmetrizing operator on all coordinates Zs^k- The wave functions 
^ = P{(j))Dxn and ^hard — ^i-P(0 — N + l)^iXn can be written in the following form, 
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n-1 N 

s=0 i<j 

N ^ 

X n , (63a) 



*ha.d = *iSYP(0 -N+1) n[ef ?*.,2 



s 

N N 



where SY is the symmetrizing operator on the electron coordinates, and 

N' 

*.,2= n (^sM-^sm)- (64) 

ki<k2 

Before doing the anti-symmetrizing or the symmetrizing in the formulas above, it appears 
that there are n different groups of electrons and there are correlations between different 
groups. The generic terms before doing the projection, for example, in the formula of ^, are 

^s,k^s,k 

(65) 



(l + ^.,fe^.,fc)^+"-' 
It will be projected to 

(0 + l)m(0 + n-l-O! 



(66) 



{<t> + n)\{l-s)\{<t>-l + s)\ii + ,^y,^^^)r 

As is a constant and is not dependent on s, Z, we can discard it in the process of the 

projection. Thus P will act in the following way (discarding constant |^^^), 

p KkF{z,,k) _ 1 1 

(1 + ^.,fe^.,ik)2+»-i (1 + Zs,kZs,k)^ ('^ - ""-^kd.Jl 

x5,%(</> + n-l-^,,fca,^_J! 

xF{zs,k). (67) 
For example, by applying this formula to = PDx2, the wave function is then given by 

AT -j^ N/2 

* = ANP n - — ^ n (</• + 1 - ^oA,.) 

i=l (1 + ZiZi) 2 fc=l 

1 N/2 N/2 

n(*^,2)' n n (^o,., - z,,,,)\ m 

S=0 fel = l fc2 = l 
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The trial wave function \l/c for filling 1 — u is related to the trial wave function \& at filling 
z/ by charge conjugation, 

fYl . . . z^+m)M^i . . . zn+m) (69) 

J i + ZN+iZN+i) 

where M is the number of particles in the state is the number of electrons in 

+ M = + 1, and ^ is the complex conjugate of Note again, If we use PDxn as the 
trial wave function the projection operator P can be dropped in Eq. (|69|). However if one 
uses \l/hard = — M + 1)'^ iXn in Eq- (|69|), then the operator P can not be drooped in 

Eq. (HD. 

One can also act D on \1/ and we will get another trial wave function of the FQH state 
at filling where z/ is the filling of the state Repeatedly acting D and C on PDxn, 
we can get the trial wave functions at all observable fillings. 

VI. THE OVERLAPS BETWEEN HIERARCHICAL WAVE FUNCTIONS AND 

CF WAVE FUNCTIONS 

We perform the calculations of the wave functions symbolically by using Maple. The 
overlaps between the hierarchical wave function and the CF wave functions are calculated. 
Some overlaps between the wave functions with or without the hard core construction are 
also calculated. The formula of the trial wave functions for the FQHE in the previous 
sections need to be normalized before we calculate the overlaps. Table || lists some overlaps 
at some fillings for a small number of electrons. E^^ means a new state formed by the 
condensation of quasielectrons of parent state and means a new state formed by the 
condensation of quasiholes of parent state In all cases listed in the table, pi is equal to 
2 for 2 > 1 in the constructions of the hierarchical states. The wave functions which involve 
D, P, C operations are the CF wave functions. 

When = 3, 4, = 6, there is only one rotational invariant state, which must also be 
the ground state. This explains why some of overlaps in the table are equal to one exactly. 
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Because of limited CPU we are allowed to use, we are only able to calculate some hier- 
archical wave functions up to 6 electrons, and some CF wave functions up to 10 electrons. 



The detailed calculations can be found in Ref. |2^. In the future, we will calculate the wave 
functions with more numbers of electrons. 

From the calculations, we conclude that, the hierarchical wave functions and the CF 
wave functions are all most the same in the case of a small number of electrons. 



VII. CONCLUSIONS 

In this paper, we present a detailed discussion about the calculation of the trial wave 
functions on the sphere. The projective coordinates are used in performing the calculations. 
A self-contained derivation of the LLs on the sphere (or any surfaces with a constant cur- 
vature) using geometrical method is also given in the paper. The many-body wave function 
in the LLL are studied and classified in the angular momentum bases. We also simplify the 
formulas for the hierarchical wave functions and the CF wave functions. 

There are many interesting things which we want to study in the future. We shall use 
theories of polynomials to study those wave functions.ii'0 It would be very interesting if 
we can obtain the polynomials explicitly for the wave functions at an arbitrary number of 
electrons. 

There is a mapping between a trial wave function in the FQHE and a wave function in an 
one-dimension space.0 Because of the existence of the mapping, one may apply the method 
used to study the Calogero model to study the trial wave functions in the FQHE, and then 
it may be possible to calculate some physical quantities from the trial wave functions at an 
arbitrary number of electrons. 
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APPENDIX: LANDAU LEVELS ON COMPACT CLOSED SURFACES 

In this appendix, we will study the LLs on general compact closed surfaces, and work 
out the LLs on the sphere as an example. 

If the magnetic field and the curvature are constant, the spectrum, the wave functions 
and the degeneracy of Landau levels (LLs) can be obtained by using a very simple geometric 



argument. cfl A self-contained presentation of the idea based on Ref. |T8| will be found in this 
appendix, and some examples will be included. 

In the case that the surface is a plane, a sphere, or a torus, the spectrum and eigen- 
functions of the LLs can be exactly solved.ii For example, the LLs on a sphere with a 
Dirac-monopole on the origin, were solved by Dirac long time ago. The problem in the case 
of the surface being an open up-half hyperbolic plane with a constant negative curvature 
was solved completely where there exist a discreet spectrum (this is the spectrum of the 
LLs) in the low-energy sector and a continue spectrum in the high-energy sector.^! 



Ref. |I8| studied the problem of the LLs on the compact closed Riemann surfaces with 
Poincare metric, and obtained the discrete low-energy eigenvalues (or LLs), their multiplicity 
and wave functions. Previous to Ref. [1^, similar problem also was studied and discrete 
low-energy eigenvalues and their multiplicity was obtained by using the results from the 
mathematical literature, for example by using Selberg trace formula (see the references 
quoted in Ref. |18[). 

Why the problem of the LLs in all those surfaces mentioned above can be solved com- 



pletely? By closely following the observation in Ref. IT3, it is quite clear that the method 



developed in Ref. IT^ can be easily generalized to the case of any constant curvature surface 
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with a constant magnetic field applied on the surface (the surface can be a compact and 
closed surface, or an open surface, for example, an open up-half hyperbolic plane), and thus 
the problem of the LLs can be solved exactly in such cases. 

We will show that, if the curvature and magnetic field are constant, we can get many 
informations about the spectrum and the degeneracy of the LLs without solving the wave 
functions of the LLs explicitly by using a simple geometric argument ,0 even though the 
surface can be a very complicated one. If the magnetic field is constant, the wave functions 
of the ground states turn out to be a holomorphic line bundle defined on the surface. If 
the curvature of the surfaces is constant too, for the high LL, the wave functions of the LLs 
are obtained by repeatedly applying covariant derivatives on some holomorphic line bundles 
(which will be specified later) . The spectrum is obtained without solving the wave functions 
explicitly and the degeneracy of the LLs can be obtained by the Riemann-Roch theorem. 
If the sections of some holomorphic line bundles can be obtained, the wave functions of the 
LLs can be obtained explicitly. 

We use two simple examples to demonstrate how to use this geometric approach to solve 
the LLs. The examples are the LLs on the sphere and the open up-half hyperbolic plane. 

1. Ground States 

We will show here that, when the magnetic field is constant, the ground states satisfies 
a first-order holomorphic (or anti-holomorphic) differential equation and the ground states 
belong to the sections of a holomorphic line bundle. 

We consider a particle on a surface interacting with a magnetic field. In complex coordi- 
nates, the metric is ds"^ = Qzsdzdz and the volume form is dv = [igzz/'^jdz Adz = g^sdx Ady. 
The natural definition of the constant magnetic field to the high genus Riemann surface is 

F = Bdv = {dzA-z - d-zAz)dz A dz, (Al) 

where B is a constant. Thus we have dzAz — dsAz = igz2B/2. If the surface is closed, 
the magnetic field is then called "monopole" field and subjected to the Dirac quantization 
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condition. The flux (0 must be an integer) is given by 27r0 = J F = BV, where V is the 
area of the surface and we assume here B > {(p > 0) for simplicity. The Hamiltonian of a 
particle on the surface is given by the following equation, 

= ^[(P. - A,){P, - A.) + (P, - A,m - A.)] 

where g'^^ = , Pz = -id^ , Pj = -id^ , = {d^ - idy)/2 , and = {d^ + idy)/2. We 
define the inner product between two wave functions as < il'i\'^2 >= / dvipi x ■02. 

Define H' = ^^(Pj — Az){Pz — Ag). H' is a positive definite hermitian operator because 
< 0|if'|'0 > > for any If H'll) = 0,0^ must satisfy (P^ — Az)^) = 0. The solutions of this 
equation are the ground states of the Hamiltonian H or H'. In the case of closed compact 
surface, the existence of the solutions of this equation is guaranteed by the Riemann-Roch 
theorem. The solutions belong to the sections of the holomorphic line bundle with the 
connection given by the gauge field. The Riemann-Roch theorem tells us that 

h%L) - h\L) = deg{L) -h + 1, (A3) 

where h is the genus of the surface, h^{L) is the dimension of the sections of the holomorphic 
line bundle or the degeneracy of the ground states of the Hamiltonian H, h^{L) is the 
dimension of the holomorphic differential {L~^ x K), where K is the canonical bundle, and 
deg{L) is the degree of the line bundle, which is equal to the first Chern number of the 
gauge field, or the magnetic flux out of the surface, 0. When deg{L) > 2h — 2, h^{L) is equal 
to zero,0 thus h^{L) = (f) — h + 1. One finds that h^{L) indeed gives the right degeneracy 
of the ground states in the case of a particle on a sphere or a torus interacting with a 
magnetic-monopole field. 

In the case of non-compact surfaces, for example an infinite plane or an up-half hyper- 
bolic plane, the flux out of the surfaces are infinite, and the degeneracy is infinite too. The 
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degeneracy of the LLs turns out be infinite. Tlius Eq. also gives correctly the degen- 
eracy, as when the flux is infinite, the equation implies that h^{L) becomes infinite. When 
the surface has a boundary, for example a disc, one would expect that Eq. ( [A3| ) is replaced 
by a new index relation given by the boundary index theory. Note that, when the flux is 
much bigger than one, the degeneracy of the ground states are approximately equal to the 
flux out of the surface. 

2. Higher Landau Levels 

We study the higher LLs in the case of the curvature of the surface being constant. When 
the curvature is constant, g^^ddlngzz = C, a Liouville-like integrable equation. For the fiat 
surface, C = 0, as in the case of a plan or a torus, the spectrum and the wave functions 
of the LLs can be completely solved. When the surface is flat, the higher LLs are obtained 
by applying successively a first order differential operator to the states in the LLL. Now we 
shall generalize such construction of the LLs in the case of a flat surface to the case of a 
curved surface. 

We consider here the closed and curved surface with constant (non-zero) curvatures. It 
is easy to generalize to the case of an open surface with a constant curvature and we will 
demonstrate it in an example in the end of the appendix. When C is not equal to zero, 
one has g^z = {^/C)ddlngz2- As the magnetic field is constant, we can fix the gauge field 
as Az = —iB'd{ln gzz)/'2, and the magnetic field F is equal to Bdv, where B = 2B'C. For 
example, in the case of the Poincare metric, ds^ = y~'^{dx'^ + dy'^) , gzs = y~'^ , and C = 1/2, 
thus B = B' . For a closed surface, by Gauss theorem, the flux out of the surface is equal 
to = B{h — l)/c = 2B'{h — 1). B' must be a rational number as is an integer. For the 
negative curvature closed surface, according to Gauss theorem, we should have h > 2. On 
the other hand, for the positive constant curvature surface, h must be equal to zero, and 
thus the surface is topologically equivalent to a sphere. Without losing any generalities, we 
assume in the following discussions that i? is a positive number. For a negative B, the wave 
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functions are the complex conjugate of the wave functions in the case of a positive B. 
For any eigenfunctions of the Hamiltonian, they satisfy 

H^Ij = EiIj. (A4) 

If the domain of z intersects non-trivially the domain of z, g^zdzdz is invariant under coor- 
dinate changes, or 

Qzzdzdz = g^gdzdz (A5) 

on the intersection of the domains of z and z. Define 

D = d-{B'/2)d\ngzz,D = d+{B'/2)d\ngzz. (A6) 

D and D are transformed as 

D = {dz/dz)U-^DU , b = {dz/dl)U-^DU (A7) 

where U{z,~z) = {%)'^' '^{%)^' ■ 

We take m = 2 in Eq. ( |A^ ) for the simphcity. The Hamiltonian can be written in the 
following form, 

H = -g''DD + {B/A). (A8) 

Thus the Hamiltonian in the domain z is transformed to the Hamiltonian in the domain z 

as 

H = U'^HU, (A9) 

and the wave function is transformed as 

tP = U-^ip. (AlO) 

Therefore tjj{dz)^' ^"^ {dz)~^' ^'^ is invariant under the transformation, and it implies that ijj is 
a differential form of type Tg,/2 ' "where we use the following notation: if F{z, z){dz) [dz] 
is invariant under the transformation, then F{z^ z) is a differential form of type . 
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The ground states are given by the solutions of the equation Dip = 0. When the curvature 
is negative, C, (p, B' are positive numbers. If > 2/;, — 2, or 5' > 1, then the number of the 
solutions is (j) — h+1 according to the previous discussions. For smaller 0, some discussions 



can be found in Ref. |T8|. In the case of compact and closed Riemann surfaces with the 
Poincare metric, the wave functions in LLL were constructed by calculating the determinant 
of holomorphic sections of some bundlcH 

When the curvature is positive, C and B' are negative numbers, and /i = as shown 
in the previous discussions. Now is equal to = —2B'. The LLL states are again given 
by the solutions of the equation Dtp = 0. As |0| > 2h — 2 = —2 {h = in this case), the 
number of the solutions is equal to|0|— /i + l = |0| + l. 

To obtain the spectrum and wave functions of the higher LLs, we introduce the covariant 
derivative,^! Vz, and its Hermitian conjugate (V^)^ = — V^, 

V, -.Tl^Ti^, ,V, = /a(7^^ (Alia) 
(V.)t: Ti ^ tU , (V.)^ = -g-'-'dg\ (Allb) 

where we call g = gzz for short. Note that D is the covariant operator Vz acting on T^/j^ , 
and D = gV^ where acts on T^ij^. The Hamiltonian can be written by using the 
covariant operators, 

H - B/A = -VzV'. (A12) 

One can verify the commutation relation, 

[V" V,]T;^ = -(m + n)C. (A13) 

Assume that ipi is a state in the higher LLs and an eigenfunction of H with eigenvalue Ei, 
then ipi = —-^'Vz'V^ipi, where ei = Ei — B /A > 0. Therefore one can write ipi = V2$(l), 
where $(1) is a differential form of type Tg//2-i- More explicitly, we have 

ipi = {d-{B'/2-l)dlng)^{l). (A14) 
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Using the relation ddlng — gC, one can show that, 

- V.V^V'i = {B' - 1)CiIj, + V,[-V,V^$(1)]. (A15) 

We first discuss the case of a negative curvature surfaces, li B' > 1, one can show that 
< '0i|V2[— VzV^$(l)] >> 0. Thus one can conclude that the states of the lowest excited 
level are obtained, if there exist $ such that V2V$(1) = 0. V2V^$(1) = leads to 

= 0. The solution of 5$(1) = is $(1) = c/-^'/2$(l) with d^{l) = 0, where $(1) 
is of the form Tb'-i- By the Riemann-Roch theorem, there exist solutions of the equation 

= for 5/ > 1, and the number of the solutions or the degeneracy of this Landau 
level is the dimension of the sections of the holomorphic bundle Tb'-i, which is equal to 
{2B' — 3) (/i — 1) if 5' > 2. The energy of this LL or the lowest excited states is 

El = ^B'C - C. (A16) 

\i B' < 1, there is only the zero'th "Landau level" or the LLL (there exists a continue 
spectrum in the high-energy sector, and the states in the continue spectrum are not called 
as the states in the LLs). 

We can generalize the above discussion to higher LLs. The wave function of the k'th LL 
is given by 

= (V,)'=$(A;) 

= {d- {B'/2 - l)dlng){d - {B' /2 - 2)dlng) 

x---{d-{B'/2-k)d\ng)^k), (A17) 

with ^{k) ^ g 2 ^[k) and d^{k) = 0. ^{k) is a differential form of the type Ts'-k- Notice 
that this construction generalizes the standard construction for the harmonic oscillator. The 
difference for the constructions of the high LLs between the case of the flat surfaces and the 
case of the curved surfaces is clear now. In the case of the surface being a plan or a torus, 
the high LLs are obtained by applying successively a first order differential operator to the 
ground states. However the situation is different when the surface is curved. for k ^ 
is not the ground state of the Hamiltonian H. 
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Using Eq. ( |A13| ), we calculate the eigenvalue of the corresponding wave function ipk , 



and it is equal to 



E. = CB'ik + l)-^^^. (A18) 



The degeneracy of the k'th LL is given by the dimension of the sections of the holomorphic 
bundle of the type Ts'-k, which is equal to {2B' — 2k — l){h — 1) when B' — k > 1. Because 
the dimension of T„ is zero when n is negative, k must not be greater than B'. Hence there 
is only a finite number of "Landau levels" . 

When B' is an integer, k can take value from to B'. When k = B', the corresponding 
$(A;) is the differential form of the type Tq. Tq is a constant function on the surface and the 
degeneracy of this LL is equal to one. For the twisted boundary conditions, which would 
physically correspond to the presence of some magnetic flux through the handles. There 
does not exist a non-zero constant function which satisfies the twisted boundary condition, 
thus the dimension of To is zero, and the degeneracy of this LL is equal to zero or there does 
not exist this LL {B' — th LL). When k = B' — 1, the degeneracy of this LL is the dimension 
of the canonical bundle Ti, which is equal to h for the non- twisted boundary condition and 
is equal to h — 1 for the twisted boundary condition (this result can be obtained by the 
Riemann-Roch theorem). B could be also an half-integer. Then k can take value from 
to B' — (1/2). When k = B' — (1/2), the degeneracy of this LL is the dimension of the 
spin bundle Ti/2- The dimension of the holomorphic sections of the spin bundle generically 
is zero for the even-spin structures and one for the odd ones (or for twisted ones). It is 
possible that B' is fractional assuming that 2B'{h — 1) is an integer, and k can take value 
from to [B'] where [B'] is the bigger integer which is smaller than B'. When k = [B'], the 
degeneracy of this LL is the dimension of bundle Tb'-[b']- B' — [B'] is a fractional number 
between and 1 and the discussions of such case can be found in Ref. |T8|. Beyond those 
LLs, little is known about the continue spectrum in the case of the complicated negative 
curvature surfaces. 

To normalize iIju, we calculate the inner product < ipkl'ipk >■ By using Eq. (|A13| ). It is 
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given by the following equation, 



<i^k\i^k> =<V^$(fc)|V^<l'(A;)> 

= < <l>(A;)|(V^)tV^$(A;) > 

k 

xl[{2B' ~k-i), (A19) 



i=l 

where the inner product < > is defined as 



< >= J dvg''^{k) x $(A;). (A20) 



The definition of the inner product between two ^{k) given in Eq. (|A2(]| ) is quite natural 

2)-k- 



because $(A;) is a differential form of the type ^(b'/2) 



If $(A;) is normalized to one, then 



^ (A21) 



is also normalized to one. 

Now we come to the case of a closed surface with a positive curvature, which is a little 
bit different from the case of a surface with a negative curvature. Now we have only h = 
according to the previous discussion. The wave function is a. differential form of type T^/j^ 
with 5' being a negative number. In the formula —VzV^-?/'i = {B'—l)C4'i+'^z[—^z'^^^{^)], 
one can show that < ipi\V zi—V zV^^) >> for any negative B'. By using Riemann-Roch 
theorem, one finds that there always exists $(1) such that V^V^$(1) = 0, which leads to 
Z)$ = 0. Therefor for any B', there exists a higher LL. One can repeat the argument to 
obtain the states in the higher LLs and obtains the full spectrum and wave functions. 

The wave functions of the states in the k'th LL are again given by Eq. ( |A17| ), with 



^(k) = g^'/^^{k) and d^{k) = 0. ^{k) is a differential form of the type T^'-k- The 
degeneracy of the k'th LL is equal to the dimension of the holomorphic line bundle Ts'-k, 
which is equal to 2{B' — k){h — 1) — h + 1 = —2{B' — k) + 1 a.s h is equal to zero. 
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The energy is again given by Eq. ( |A18|) . However, a higher LL has a higher degeneracy 
and the number of the LLs is infinite in such case. Instead, in the case of a surface with 
a negative curvature, a higher LL has a smaller degeneracy, and the number of the LLs 
is finite. From Eq. ( |A18| ), one notices that, in the case of a positive curvature surface, the 
energy gap in the neighboring LLs increases when the level increases, and in the case of a 
negative curvature surface, the energy gap in the neighboring LLs decreases when the level 
increases. 

It is easy to generalize the above discussions to non-compact surfaces, and we will work 
out an example in the following discussion. 



3. Examples 

a. Upper half hyperbolic surface 

We consider that the surface is a upper half hyperbolic surface (also see Comtet and 
Dunne in Ref. ^9]). In the projective coordinates, the metric g is written as (jz^ip-, where 
1 2; I < 1 , The other quantities are, C = 2,B = 4B' ,andA, = -ijff^- The wave functions 
are given by Eq. ( |A17| ) . As the wave functions of the LLs shall be normalizable (opposite to 
the wave function of a state inside the continue spectrum), < ipkl'^^k > shall be normalizable. 
A normalizable < ipk\4'k > is equivalent to a normalizable < $(/c)|$(A;) >. A normalizable 
< $(A;)|$(A;) > leads the condition B' — (1/2) > A; > 0. ^{k) is given by function 
where I is a non-negative integer. Thus the degeneracy is infinite for every LL. This is 
consistent with the Riemann-Roch theorem as the flux out of the surface is infinite. Finally, 
the energy is given by Eq. ( [A18| ). 



b. Sphere 

Another example is that the surface is a sphere. In the projective coordinates, the 
metric g is written as g = • The other quantities are, = ij]^-^ , C = —2 , and 
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B = 2B'C = —4B'. Thus the flux = —2B' is a non- negative integer according to the Dirac 
quantization condition (note we always assume i? > in this paper). The wave functions 
are again given by Eq. ( |A17[ ), the energies are given by Eq. ( [A18| ), and $(fc) is given by 
g't>/4:^i_ The normahzable condition leads to I = 0,1, ■ ■ ■ ,2k + (p. Thus the degeneracy of the 
k'th LL is equal to 2A; + + 1. The degeneracy can be also obtained by the Riemann-Roch 
theorem and the result is consistent with the result obtained by requiring the wave functions 
being normalizable. In this way, we obtain the full spectrum and all wave functions on the 
sphere. 

The wave functions at the n'th Landau level (n = is the lowest Landau level) are given 



From previous discussions, we can easily flnd the inner product < ipnAi^n^i > is equal to 



However inside the paper. The deflnition of the inner product is different from the defi- 
nition in the appendix. The inner product in the paper is defined as 



by Eq. (g). 



TT 



;!(0+2n-O! 



. The inner product is as previously defined, < '?/'i|'?/'2 >= / dvipi x -02, where 



{<f>+2n+l)((f>+nY. 





(A22) 



As dzdz = 2dxdy, thus < 'ipn,i\'ipn,i > is given by the following formula. 



< i^n,l\i^n,l >= 27r 



/!(0 + 2n-/)! 



(A23) 



(0 + 2n + l)(0 + n)!' 



This formula is used in the paper. 
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TABLES 

TABLE I. In this table, we list the number of rotational invariant states at various v and a 

small number of electrons. Ni is the dimension of the total Hilbert space (in the LLL) and A^ is 
the number of the rotational invariant states. 
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N 


At 


A^ 


^ (formula) 


4> 

T 


N<t> 
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2 
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5 
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|A-4 
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12 
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3 




11 


33 
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8 




16 


64 


2 
7 


4 


43 
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^A-2 


12 


24 
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1.242 
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19 


57 
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46.029 
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26 
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2 
9 


4 


43 
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|A-6 
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12 


24 
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13 




21 


63 
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139.143 
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^■N - 4 
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18 


36 
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1.229.093 
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4 
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i^A-8 
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18 


36 
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17.002 
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31 


93 
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1.137 




44 


176 
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8 
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16 
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3 
11 
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2.137 


13 


^A-1 


21 


63 
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610.358 


506 




32 


144 


3 
17 
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17.002 


34 


fA-3 


31 


93 


5 

17 


4 


33 


2 


5 5 


11 


22 
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184.717 


217 




28 


126 
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TABLE II. The overlaps between the hierarchical wave functions and the CF wave functions 

at some lilliugts for a small numV)(;r of elcclrons. 
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overlap 


2 
5 




< PDx2\E^3 > 


< *iP*iX2|£^*3 > 


< PDx2\'^iE'^i > 


< *iP*iX2|*l£^*l > 
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.9993234149 
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.9999456457 
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1 


1 


1 


1 



v 


N 
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v 


N 
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2 
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< E'^sl'^iE^l > 


< PL>X2|^'l-P^lX2 > 


2 
7 




< H^3\DC^3 > 
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.9996479001 


.9999288987 
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.9993762574 
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1 
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2 
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< E^5\DPDx2 > 


< E^5\PD'^X2 > 


< DPDx2\PD^X2 > 
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.9999614869 


.9999614869 


1 
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3 
11 




< EH^3\DCPDx2 > 
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< H^5\D'^C^3 > 
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.9996522383 
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< E^i\DPD'^x2 > 
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< HH-^alDCDC^a > 
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.9999999997 
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